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Abstract. The eigenstates of d-dimensional quasicrystalhne models with a separable 
Hamiltonian are studied within the tight-binding model. The approach is based on mathematical 
sequences, constructed by an inflation rule V = {w s,s ^ sws^~^} describing the 
weak/strong couplings of atoms in a quasiperiodic chain. Higher-dimensional quasiperiodic 
tilings are constructed as a direct product of these chains and their eigenstates can be directly 
calculated by multiplying the energies E or wave functions ^' of the chain, respectively. 

Applying this construction rule, the grid in d dimensions splits into 2^*"^ different tilings, 
for which we investigated the characteristics of the wave functions. For the standard two- 
dimensional labyrinth tiling constructed from the octonacci sequence (6 = 2) the lattice breaks 
up into two identical lattices, which consequently yield the same eigenstates. While this is not 
the case for 6 7^ 2, our numerical results show that the wave functions of the different grids 
become increasingly similar for large system sizes. This can be explained by the fact that the 
structure of the 2''"^ grids mainly differs at the boundaries and thus for large systems the 
eigenstates approach each other. This property allows us to analytically derive properties of the 
higher-dimensional generalized labyrinth tilings from the one-dimensional results. In particular 
participation numbers and corresponding scaling exponents have been determined. 



1. Introduction 

Since the discovery of quasicrystals research aims at understanding the physical properties of 
these materials. Today many exact results are known for one-dimensional quasicrystals [HE] but 
the characteristics of wave functions in two or three dimensions have been clarified to much lesser 
degree. In particular, they have mainly been studied by numerical calculations of finite systems 
or periodic approximants due to the missing translational symmetry of these materials. For 
instance numerical results for the Penrose tiling, the Ammann-Kramer-Neri tiling or octagonal 
tiling models have been determined but for relatively small systems making it difficult to estimate 
the behavior for macroscopic system sizes j3ll3]. Further, it is possible to calculate some exact 
eigenstates of the tight-binding Hamiltonian on the Penrose tiling [5l[6], which is a first step to 
understand the nature of eigenstates and properties of quasicrystals in higher dimensions. 

In order to get a better insight into the properties of wave functions in quasicrystals, 
another approach employs models of d-dimensional quasicrystals, which can be built from 
d separate one-dimensional systems yielding straightforward solutions. This allows us to 
study very large systems in higher dimensions based on the solutions in one dimension 
[7]. Unfortunately, these structures do not occur in real quasicrystals, however, they can 
be artificially constructed. Examples are conducting nanowires, coupled nanomechanical 



resonators, or photonic quasicrystals [TJ [H]. The best studied model within this approach is 
based on the Fibonacci sequence corresponding to the golden mean [UE], whereas in this paper 
we concentrate on the other metallic mean quasiperiodic sequences UHl E] , which have been 
less investigated so far. 

At first we introduce the quasicrystal model based on metallic mean sequences and discuss 
the symmetry properties of the eigenstates in Sec. [2] This approach is expanded to higher 
dimensions for the labyrinth tiling in Sec.jsj In general, the grid splits into 2*^"^ different tilings 
in d dimensions but the results indicate that the wave functions of these tilings approach each 
other for large systems. We give a reasonable argument for this behavior as well. This property 
can be used to calculate or even analytically derive characteristics for generalized labyrinth 
tilings. As an example, Sec. |4] comprises results for the participation ratio of the eigenstates and 
their scaling behavior in one, two, and three dimensions. The paper is concluded in Sec. [5] 

2. Eigenstates and wave functions of metallic mean chains 

The metallic mean sequences for a parameter b are defined by the inflation rule V = {w — )• 
s, s — )• sws^^^} based on an alphabet containing the symbols {s, w} with the starting symbol w. 
After a iterations we obtain the ath order approximant Ca of the quasiperiodic chain. Regarding 
the length fa of this approximant Ca we can also find a recursive rule fa = hfa-i + fa-2 with 
/o = /i = 1- The ratio of the lengths of two successive iterants in the limit lima_^oo fa/fa-\ = A 
is given by different metallic means depending on b with a continued fraction representation 
A = [6] = [6, 6, 6, ...]. This leads to the well known Fibonacci sequence for 6 = 1 with the golden 
mean Aau = [!] = (! + \/5)/2, while 6 = 2 results in the octonacci sequence with silver mean 
AAg = [2] = 1 + \/2 and 6 = 3 corresponds to the bronze mean Abz = [3] = (3 + \/T3)/2 pT] . 
For the orthogonal basis states |/) the tight-binding Hamiltonian can be written as 

fa fa + l 

H = Y,\i)%i+^{i + i\ + Y,\i)^dA ■ (1) 

1=1 1 = 1 

This can be interpreted as an electron on a discrete grid which can hop from one vertex of 
this graph to any of its neighboring vertices and the aperiodicity of the grid is given by the 
quasiperiodic sequence with the symbols w and s corresponding to weak and strong bonds. 
Hence, the hopping parameter t is determined according to the sequence Ca with = 1 for a 
strong bond and tyj = v for a weak bond (0 < f < 1). Further, usually one assumes zero on-site 
potentials (e/ = 0) for these models so that no vertex is energetically preferred [21 [10] . 
The discrete energy values E"^ and corresponding wave functions |^'^) = E/=t^*ilO 

are 

solutions of the time-independent Schrodinger equation for the quasiperiodic systems 

^ I**) = E' \¥) =^ E'¥ = ti^i,i¥i_, + tM+i^«+i • (2) 

The results show that the eigenvalues are symmetric with respect to 0. However, we have 
to distinguish between even and odd system sizes Na = fa + ^- for even Na all energy values 
E have a symmetric value —E but for odd Na there is one state E^'^ = 0, which has no 
corresponding state. Additionally, also the eigenfunctions show a symmetry. The eigenstates ^ 
with an eigenvalue E and ^I' with the opposite eigenvalue —E only differ by an alternating sign 
depending on the position / with ^ 

^« = i-iy^i . (3) 

Further, for odd Na = 2M — 1 the eigenvector associated to the eigenvalue E^^ = has a 
special structure, namely vanishes either on all odd or on all even sites I (cp. [lO]), i.e., 

^fmod 2.0 = or *fCod2.1 = 0. (4) 



3. Eigenstates and wave functions of the generalized labyrinth tiling 

The generalized labyrinth tiling is based on two quasiperiodic chains perpendicular to each 
other. Hence, the labyrinth tiling can be expressed as the Euclidean product £a = Ca x Ca of 
the two chains. The bonds are given by the diagonals of the constructed grid and depending 
on the starting point the grid decomposes into two separate ones, cp. Fig. [T]^a). Connecting all 
vertices with even index sum, we obtain the finite ath order approximant Ca of the generalized 
labyrinth tiling C. Analogously, all vertices with odd index sum connected by diagonal bonds 
also form a labyrinth tiling £* [9l [10]. Both tilings are shown for different inflation rules in 
Fig. [TJ For the octonacci sequence they are identical if Ca is mirrored at the centre of the x- or 
y-axis or rotated by 90°, i.e. i2* is dual to Ca- For other inflation rules Ca and £* are slightly 
shifted against each other and mainly differ at the boundaries as shown in Fig. [l]^b) and[T]^c). 
These similarities originate from hidden mirror symmetries in the quasiperiodic chains. While 
the octonacci chain is perfectly mirror symmetric, for the other chains mirror symmetry can be 
achieved by neglecting a few symbols at one end of the sequence (e.g. for the Fibonacci chain 
the last two symbols on the right) and by interchanging two consecutive bonds for the case 6 > 2 
(i.e. a phason flip). 

The eigenstates of the generalized labyrinth tiling are constructed as the product of the 
eigenstates of two one-dimensional chains, i.e. E^^ = E^^E"^^ and ~ ^'p^m [IQj- The indices 
i and j enumerate the eigenvalues E in ascending order and the corresponding wave functions 
^I', and / and m represent the coordinates of the vertices. However, due to the symmetries of 
the eigenfunctions of Eq. ([s]) , some of the wave functions and related eigenvalues are 
identical. Then only one of them is allowed to be considered. In two dimensions all combinations 

{E'^\l<i<Na/2M<J<Na] (5) 

are valid for even chain lengths Na- Further, the wave functions have to be normalized. For 
odd Na the selection of the eigenstates is more complicated due to the special structure of the 
eigenfunction for E^^ = 0. This results in the two-dimensional energy values 

{E'^ \ {I < i < M M < j < Na) y {i = M M < j < M)] . (6) 

In three dimensions we follow the same considerations as in two dimensions. The grid now 
decomposes into four different grids depending on the choice of the starting point because the 




(a) silver mean model C^^* (b) golden mean model C^^ (c) golden mean model C^' 



Figure 1. Labyrinth tilings Ca and for different inflation rules: (a) visualizes the construction 
of Ca and £* for the octonacci chain; (b) shows C^^ and (c) for the golden mean model, 

where Ca is mirrored at the centre of the y-axis to show the similarities of both grids. 



bonds correspond to one of the four main diagonals of a cuboid. Without loss of generality we 
consider only the labyrinth tiling C'^ = Ca x Ca x Ca, for which all vertices have either all even 
or all odd indices /, m, and n. 

We analogously compute the eigenvalues and vectors in three dimensions. Here they are based 
on three independent quasiperiodic chains, where each one is perpendicular to the other two 
chains. The wave functions are now given as a product of three one-dimensional eigenfunctions 
with ^^i^^ ~ ^Y'^^^^^ , where i, j, and k label the energy values and m, and n the vertices 
of the tiling. For even Na the eigenvalues E'^^^ = E^E^E^ of the generalized labyrinth tiling C^'^ 
in three dimensions are given by the set 

{e'^'' \l<i< Na/2 A 1 < J < Na/2 A 1 < A; < iV,} . (7) 

This equation results from the symmetry of the wave functions from Eq. ([s]) . Due to the structure 
of the eigenfunctions to the eigenvalue i?*^ = 0, for odd Na the allowed combinations are given 
by 

{^iifc I {1 < i < M Al < j < M Al < k < Na)y 

{i = M Al < i < M Al <k < M)\J {I <i < M A j = M Al <k < M)} . (8) 

The characteristics of the wave functions <^ for the different tilings are interesting. Figure [2] 
shows several wave functions of the generalized labyrinth tilings Ca and in two dimensions. 
For the silver mean model the wave functions are identical for Ca and £* if one of them is 
mirrored at the centre of the x- or y-axis, cp. Fig. ^a) andpfd). Both wave functions are also 
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Figure 2. Similarities of wave functions for the two different generalized labyrinth tilings Ca 
and C* shown for the silver mean model and for two different approximants of the golden mean 
model. The point size and color encode the probability |$p with a logarithmic scale, cp. |6]. 



identical if rotated by 90° taking into account that the indices i and j of one of the eigenstates 
have to be interchanged. For other inflation rules the wave functions are not necessarily identical 
for small approximants due to the different boundaries of the lattice. Nevertheless, they are quite 
similarly distributed over the lattice as e.g. displayed in Fig.[2]^b) and[2]^e) and several of the wave 
functions show even a symmetry with respect to the x- or y-axis. For larger approximants the 
wave functions approach each other because the influence of the boundaries and of the phason 
flip for b > 2 vanish, see Fig. and^f). This characteristic can be helpful for analytical 
considerations for properties of the general labyrinth tiling. 



4. Scaling behavior of participation ratios 

The participation ratio p reflects the fraction of the total number of sites for which the probability 
measure of the wave function |$rP is significantly different from zero and, hence, provides 
information about the degree of localization of the wave function <I>r- For the discrete positions 
r (all valid combinations of the indices /, m, and n) it is defined via the inverse participation 
number P and the number of states Va (= N^/2'^~^ for even Na) in d dimensions as 



P($) 
Va 



1 



(9) 



The participation ratio Va) of a state scales with a power law in the number of sites 
as p{^,Va) ~ Va where the exponent 7 is connected to the nature of the eigenstates [3j. 
The wave function of a localized state is characterized by 7 = 1, and 7 = corresponds to an 
extended state. For intermediate values of 7 (0 < 7 < 1) we obtain fractal eigenstates, which 
are neither extended over the whole system nor completely localized at a certain position and 
show self-similar patterns [6l [121 US] • 

Usually, the participation ratios in a small energy interval are distributed over a certain range 
resulting in a snow-flurry-like appearance when plotted. However, the average participation ratio 
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Figure 3. Scaling behavior of (p) for the silver mean model in d = 3. The dots indicate {p) 
for different system sizes and the lines correspond to a mean least squares fit, where the slope 
equals the scaling exponent 7. 
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Figure 4. Scaling exponents 7 of the average participation ratio {p) over the complete energy 
spectrum in one, two, and three dimensions for different metallic mean models. 

is almost identical for states at the centre and the edge of the energy spectrum for the present 
model |12] . For other quasiperiodic tilings differences can be observed, e.g. the participation 
ratios of eigenstates of the Penrose tiling, the generalized Rauzy tiling and octogonal tiling 
models are smaller at the centre than at the edges of the spectrum, whereas e.g. the Ammann- 
Kramer-Neri tiling shows an inverse behavior [3l HI [Tl| . In order to compare the characteristics 
of the eigenstates for different dimensions we compute the scaling behavior for the average 
participation ratio {p{^,Va)) = ^J2sPi^^^^a)- Results for the silver mean model for three 
dimensions are shown in Fig. [3j The average participation ratio scales very well with the system 
size. In Fig.|4]we show the dependency of the scaling exponent 7 on the coupling constant v for 
different metallic means in one, two, and three dimensions. The plot indicates that the scaling 
exponent 7 is independent of the dimension, which can be also analytically proven for the limit 
of infinite systems. The deviations for the different dimensions, especially for d = 3, for the 
golden and bronze mean model are caused by the higher computational efforts allowing only the 
consideration of comparatively small systems. For the octonacci chain participation ratios of 
higher-dimensional tilings are given as a product of the participation ratios of the chains |15j . 

The similarities of the wave functions of the different general labyrinth tilings allow us to 
find good approximations for the calculation of certain properties such as participation ratios. 
In Fig. [5] we compare the approximate results obtained by the product approach with the exact 
results and find that the deviations become very small already for system sizes with Na ~ 20. 

5. Conclusion 

In this paper we investigated the eigenstates for the generalized labyrinth tilings. We presented 
the construction rules for the eigenstates in two and three dimensions and found that for inflation 
rules with b 2 the wave functions of the 2'^~^ different tilings in d dimensions approach each 
other with increasing system sizes in analogy to the occurrence of identical wave functions for 
6 = 2 for the different grids. This property is useful for approximate calculations of certain 
properties as shown for the participation ratios but also for analytical considerations. 
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(a) golden mean model {b = 1) (b) bronze mean model (6 — 3) 

Figure 5. Difference of the average participation ratio for the exact numerical result with the 
results obtained using a product approach for different generalized labyrinth tilings in three 
dimensions. 
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